GRB Explosion in a Galactic Supershell : GRB 971214 
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Among a number of gamma-ray bursts whose host galaxies are known, GRB971214 
stands out for its high redshift(z = 3.25) and the Lya emission line having a P- 
Cygni type profile, which is interpreted to be a direct consequence of the expanding 
supershell. From a simple Voigt-fitting analysis, we estimate the expansion velocity 
of the supershell and the neutral column density. The theory on the evolution of 
supernova remnants is used to propose that the supershell is at the adiabatic phase, 
with its radius R = 10-Egg 2 pc, its age t = 4.7 X 10 3 E^ 2 yrs, and the density of 

the ambient medium n\ = 5.4£ 53 1//2 cm — 3 , where E53 = E/10 53 ergs. I estimate 
the kinetic energy of the supershell to be E^ = 7.3 X 10 52 -E53 ergs. I provide an 
observational evidence that the gamma-ray burst can occur in a giant H II region 
whose environment is similar to that in star-forming galaxies. 



1 Motivation 

Gamma ray bursts (hereafter GRBs), located at cosmological distances, form 
a group of the most luminous objects in the Universe. A number of GRBs 
were observed with their host galaxies, of which intensive spectroscopic and 
imaging observations have been performed using the Hubble Space Telescope 
and the Keck telescopes. 

GRB971214 is one of such objects with a very high redshift z > 3 and also 
is worth a particular attention in following points. 

Firstly, even after the optical transient region had faded away, we can 
marginally detect a bright spot in the HST image, and the continuum and 
emission line fluxes from this spot overwhelm those from the remaining part of 
the host galaxy. Secondly, the ultraviolet spectrum of GRB971214 illustrated 
in Figure 7.1 shows a flat UV continuum which is often found in star forming 
galaxies. In addition, the Lya emission line has a black absorption trough in 
the red part of the emission peak. These facts imply that the UV spectrum is 
formed in a star forming region which is surrounded by a thick and expanding 
medium of neutral hydrogens. We note that the location of the GRB afterglow 
coincides with the star forming region or the bright spot, which leads to the 
proposal that the GRB occurs in a star forming region, and in this Letter we 
deduce the physical environment of the GRB from its spectrum. 
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The P-Cygni type Lya emission in the spectrum of primeval galaxies has 
been often attributed to an absorption effect by a galaxy not associated with 
the primeval galaxy but intervening accidentally in the line of sight. It has 
been regarded as a damped Lya absorption that occurs in the vicinity of the 
source galaxy. In order to check this possibility, we calculate the probability 
for observing an intervening galaxy in front of the GRB host galaxy, which is 
none other than the optical depth for seeing a galaxy between the GRB host 
galaxy at z = 3.425 and the place that corresponds to t> cxp in the Hubblc's 
expansion law. The optical depth is simply expressed by 



where the comoving volume number density of normal galaxies at z = 0, n g ~ 
0.02/i 3 Mpc~ 3 (Im 1995), and the path length L is estimated to be L — v cxp /H 
with H being the Hubble constant at the redshift and given by 



where the cosmological density parameters are SIm = SttGpo/3Hq, tt\ = 
A/3Hq, and the Hubble parameter H = 65km s _1 Mpc -1 . Here, the cross 
section is given by a — 7r(riol0/i _1 kpc) 2 , where rio is the typical galaxy size 
in units of lOkpc. A direct substitution yields the optical depth r = 0.0023r^ 
for M = 1/3 and Q A = 2/3, and r = 0.0025r? for Q M = 0.3 and fi A = 0. 
This indicates that if the average size of galaxies at z ~ 3 is not large, then 
it is highly improbable that the damped absorption in the spectrum of host 
galaxy of GRB971214 is formed by a galaxy intervening accidentally. 

This leaves us to consider an alternative hypothesis, according to which 
the P-Cygni type profile of Lya is formed by the the expanding supershell 
that surrounds the star forming region in GRB971214 and is the remnant of 
the GRB precedent to GRB971214. In order to check this possibility, we now 
calculate a number of GRB event in a galaxy at z ~ 3.4. Assuming the super- 
nova rate is proportional to the star-forming rate, Sadat et al. (1998) calcu- 
lated the supernova rate at z « 3.4, Y S n ^ 0.011 x 10 6 h% 5 SNeMyr~ 1 Mpc -3 , 
where Hq = 65 h^km s _1 Mpc 1 . Accepting the concept that the GRB rate 
traces the massive star formation rate, Woods & Loeb (1998) showed T GRB ~ 
10" 6 r SAr , where il M = 0.3, n A = 0.7, and H = 65 h 65 km s _1 Mpc _1 . There- 
fore, Tgrb(z ~ 3.4) ~ 0.011 ft,g 5 Myr" 1 Mpc -3 . Adopting the number density 
of galaxies at 2.0 < z < 3.5, $* = 1.76 x 10~ 3 /i| 5 Mpc" 3 (Pozzetti et al. 1998), 
we can get the GRB rate per a galaxy at z s=s 3.4, Tcrb — 6Myr~ 1 . Therefore, 
the number of GRB events per a galaxy during 10 4 yrs is Nqrb = 0.06. More- 
over, the beaming factor, if exists, can increase the event rate by another factor 



t = n g (l + zfaL, 



(1) 



h = H [n M (i + z) 3 + n A }^ 2 



(2) 
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of ten to hundred and the number of GRB events per a galaxy during 10 4 yrs 
is Nqrb ~ 0.6 — 6, which makes our supershell hypothesis more probable and 
alternative suggestion. 

In this talk, we adopt the hypothesis of GRB-driven supershell and per- 
form a profile fitting analysis to derive physical parameters characterizing the 
expanding supershell of neutral hydrogens surrounding the star forming region 
of the GRB host galaxy. 

2 Images and Spectrum of GRB971214 

GRB971214 was detected at 9 UT, December 14, 1997 (Heise et al. 1997), and 
its optical counterpart twelve hours after the burst (Halpern et al. 1998). With 
a total fluence of 1.09 x 10 -5 ergs cm~ 2 (Kippcn et al. 1997) and the measured 
rcdshift of z = 3.418 (Kulkarni et al. 1998), its energy release is estimated to 
be ~ 3 x 10 53 ergs in 7-rays alone under the assumption of isotropic emission, 
Oo = 0.3, and Ho = 65km s _1 Mpc . 

In Figure 1 is shown its ultraviolet spectrum rcdshiftcd to the optical 
band and obtained by Kulkarni et al. (1998) using the Keck telescope. It is 
characterized by a flat UV continuum that is typically found in the spectra 
of star forming regions. It is also seen that the Lya emission has a P-Cygni 
type profile, which is frequently observed in the astronomical objects near 
or far (Lee & Ahn 1998). Lee and Ahn (1998) proposed that the P-Cygni 
type Lya line is formed when the Lya photons emitted in the central super 
star-cluster are radiatively transferred in a HI supershell that are optically 
thick and expanding. Hence, we assume that the photon source is the H 
II region that may contain 10 4 O stars. According to Marlowe et al.(1995), 
nearby starbursting dwarfs are inferred to contain a similar number of OB stars, 
when considering Ha luminosity. So this can be thought to be neither entirely 
new nor extreme assumption for galaxies of higher rcdshifts. Furthermore, it 
appears less plausible that the Lya emission arises from the medium ionized 
by shocks produced by supernovae or hypernovae. This is because the number 
density of the inner region is not sufficiently high to give the recombination 
time scale < 10 5 years. 

3 Is the Photon Source Surrounded by the GRB Remnant? 

3.1 Interpretation of the P-Cygni Absorption 

In this work we will consider the shell hypothesis, and derive the physical 
properties of the expanding neutral medium from the observed Lya absorption. 
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We assume a Gaussian profile for the unobscured Lya emission and con- 
volve it with a Voigt function with the center displaced by the expanding 
velocity that will be determined by the fitting procedure. In principle, the 
effect of the frequency redistribution by back-scatterings should be considered. 
However, we neglect this effect in this paper, because the S/N ratio and the 
resolution of the spectrum are not sufficiently good. For the continuum level, 
the blue part of Lya, which is more prone to extinction, is extrapolated from 
the red portion of the spectrum given by Kulkarni et al. (1998). They quote 
F v = 17A{v/vn) a nJy with a = —0.7 ±0.2, where F v is the spectral density at 
frequency v and vr = 4.7 x 10 14 Hz, the central frequency of the R band. 

We show the result in Figure 1, where the dotted line represents the best 
fit profile, the solid line the observed profile, and the horizontal solid line the 
continuum level. The best fit expansion velocity of the supershell relative to the 
H II region is determined to be w OX p = 1500 km s~\ and the best fit line center 
optical depth to = 6 x 10 6 , which corresponds to Nm = 10 20 cm -2 . The best fit 
Lya profile has the width of a — 5 A and the line center flux /(A = 5280.8 A) = 
0.675 ^iJy, which gives the unobscured flux to be 9.1 x 10 -18 erg cm~ 2 s^ 1 and 
the systemic redshift z = 3.425. 

Assuming a standard Friedman cosmology with Hq — 65 km s _1 Mpc -1 
and fio = 0.3, the luminosity distance dh = 9.7 x 10 28 cm. Considering 
the Galactic extinction, the unobscured Lya flux is corrected to be Fz, ya = 
(1.5±0.7) x 10~ 17 crgcm~ 2 s _1 , where the observational error given by Kulka- 
rni ct al. (1998) is introduced. Therefore, for the assumed cosmology, the 
Lya line luminosity L^ ya = (1.8 ± 0.8) x 10 42 crg s _1 . If there is no internal 
extinction in the interior of the Lya source, this corresponds to n e = (1.4 ± 
^(Tt) - 5 ^)- 1 - 5 cm" 3 or n e = (40 ± lO)^)^^)"^ cm"*, of 
which the ionization can be maintained by ~ 10 4 05 stars as the ionizing 
source. 

From the Lya luminosity, we can estimate the star-formation rate (Thomp- 
son, Djorgovski, & Trauger 1995) to be R SF = (7 ± 3)M Q yr _1 , with both 
the internal and the Galactic extinction being corrected. This is consistent 
with the star forming rate given by Kulkarni et al. (1998) as a lower limit, 
Rsf = 5.2 M Q yr _1 which was obtained from the rest-frame continuum lumi- 
nosity at 1, 500 A. 

3.2 Physical Configuration of the Supershell 

We consider the dynamical evolutionary model of the supernova remnant to 
derive the physical quantities of the shell. According to Woltjer (1972), the 
supernova remnant has four evolutionary phases, that is, the free expansion 
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phase, the Sedov- Taylor or adiabatic phase, the snowplow or radiative phase, 
and finally the merging or dissipation phase (see also Reynolds 1988). 

According to Woltjer(1972), the radiative phase begins roughly when the 
expansion velocity of the shell becomes 

" = 300 (^V^V-* Y'" (3) 
Vlcm- 3 / V 10 er S s / 

where n\ is the number density of the ambient medium and E is the initial 
explosion energy. Since the expansion velocity v — v cxp = 1500 km s _1 of the 
supershell exceeds the velocity in the radiative phase by a large margin, we 
propose that the supershell is in the adiabatic phase, which is described by the 
Sedov solution. 

According to the Sedov solution in a uniform medium of number density 
m in which we have the relation m = 3N/R, 

/ F \ 1/4 

R = 0.92 — — t 1 ' 2 cm, (4) 
\m H N J 

/ E \ 1/4 

^°- 37 (^) *- V2 -s-\ (5) 
»i=3.2 (^) 1/ \-/ 2cm -3 j (6) 

where v is the expansion velocity of the supershell, R the size of the supershell, 
E the initial explosion energy, N the column density of the supershell, m# the 
hydrogen mass, and t the age of the shell. 

Using the values v — v CKp — 1500 km s -1 and N = 10 20 cm~ 2 , we get 

t = 4.7xl0 3 (|g) 1/2 yrs, (7) 
/ F \ 1/2 

m = 5.4 g) 1/2 cm" 3 , (9) 

where N 20 = N/10 20 cm- 2 and E 53 = E/10 53 ergs. 

From these values we estimate the total kinetic energy of the expanding 
supershell given by E k — 2irR 2 N m 'mHV 2 = 7.3 x 10 52 E 53 ergs. It is also 
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noticeable that this kinetic energy can be comparable to those of the galactic 
supershcll including those in Our Galaxy, NGC 4631, and M101 (Heiles 1979, 
Rand & van der Hulst 1993, Wang 1999). Furthermore, the P-Cygni Lya lines 
of the primeval galaxies also show the similar energy scale. Thus, we propose 
that the supershcll can be the remnant of a hypernova or a GRB that had 
exploded earlier than GRB971214. 

4 Summary and Implications 

We have studied on the formation of P-Cygni type Lya in the spectrum of the 
host-galaxy of GRB971214, and found that there are at least three components 
in the system, i.e. the parsec scale remnant of GRB971214 itself, a giant H II 
region, and a supershcll surrounding it. 

The giant H II region from which Lya emission originates is photoionized 
by a super stellar cluster whose total ionizing photons correspond to those emit- 
ted from about 10 4 05 stars. The existence of the P-Cygni type absorption 
plausibly implies that there exists a supershcll surrounding the H II region. By 
a profile fitting procedure, we found out that the shell is expanding with a ve- 
locity of w xp = 1500 km s _1 and its neutral column density iVm = 10 20 cm~ 2 . 
We also revised the redshift of the Lya emission source to be z = 3.425, and 
its unobscured Lya luminosity to be Ll V u = (1-8 ± 0.8) x 10 42 erg s _1 , which 
gives a more reasonable star formation rate to be Rsf = (7 ± 3) M Q yr _1 . 

We also applied the theory on the hydrodynamical evolution of supernova 
remnants to the supershell surrounding GRB971214. Assuming a reasonable 
scale of the initial explosion energy of the supershcll, we propose that the 

1 /2 

supershell is at the adiabatic phase, with its radius R = 18 pc, its age 
t = 4.7 x 10 3 E^ 2 yrs, and the number density of the ambient medium n\ — 
5.4 £? 53 1//2 cm~ 3 , where E53 — E/10 53 ergs. And we estimate the kinetic energy 
of the supershell to be Ek = 7.3 x 10 52 i?53 ergs. 

Recently one of the most debated suggestions is the hypernova conjecture, 
according to which gamma-ray bursts occur in star forming regions. Our results 
strongly favor this model and more concrete evidence is expected as the sample 
of GRB spectra showing Lya emission becomes statistically significant. If 
gamma-ray bursts occur in star-forming region, then we can expect the density 
of ambient media is very large and exceeds 10 — 100 cm~ 3 . However, recent 
multi- wavelength analysis of selected gamma-ray burst afterglow by Panaitescu 
& Kumar (2001) have shown that the ambient density is as low as 10 _1 to 
10 -3 cm -3 , and perhaps even less. Scalo & Wheeler (2001) argued that low 
ambient density regions naturally exist in areas of active star formation as the 
interior of superbubbles. This idea had been already presented in Ahn (1998), 
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where the observational evidence of pre-existing superbubble in the host galaxy 
of GRB971214. 
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Figure 1: Profile of the Lya emission line and its P-Cygni type absorption (solid line) with 
the best fit profile (dotted line). The emission has a Gaussian profile whose width is 5.0 A, 
line center flux /(A = 5380.8 A) = 0.675 /iJy. The line center frequency A = 5380.8 A 
gives z = 3.425. The P-Cygni type absorption is best fitted by the column density of the 
supershell iVni = 10 20 cm — 2 and the expansion velocity « oxp = 1500 km s — 1 . The dotted 
vertical line denoted by 'A' is the solid line center of our result, the line denoted by B is the 
location of the P-Cygni absorption trough, and the 'C vertical line is provided by Kulkarni 
et al. (1998). Comparing with Lye? we also show SiIIA1190 line, but the S/N ratio is not 
good. 
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